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We formed the GEnetics of Nephropathy-an International Effort 
(GENIE) consortium to examine previously reported genetic asso- 
ciations with diabetic nephropathy (DN) in type 1 diabetes. GENIE 
consists of 6,366 similarly ascertained participants of European 
ancestry with type 1 diabetes, with and without DN, from the 
All Ireland- Warren 3-Genetics of Kidneys in Diabetes U.K. and Re- 
public of Ireland (U.K. -R.O.I.) collection and the Finnish Diabetic 
Nephropathy Study (FinnDiane), combined with reanalyzed data 
from the Genetics of Kidneys in Diabetes U.S. Study (U.S. 
GoKinD). We found little evidence for the association of the 
EPO promoter polymorphism, rsl61740, with the combined phe- 
notype of proliferative retinopathy and end-stage renal disease in 
U.K.-R.O.I. (odds ratio [OR] 1.14, P = 0.19) or FinnDiane (OR 1.06, 
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P = 0.60). However, a fixed-effects meta-analysis that included the 
previously reported cohorts retained a genome-wide significant 
eissociation with that phenotype (OR 1.31, P = 2 X 10"^). An 
expanded investigation of the ELMOl locus and genetic regions 
reported to be associated with DN in the U.S. GoKinD yielded only 
nominal statistical significance for these loci. Finally, top candi- 
dates identified in a recent meta-analysis failed to reach genome- 
wide significance. In conclusion, we were unable to replicate most 
of the previously reported genetic associations for DN, and signif- 
icance for the EPO promoter association was attenuated. 
Diabetes 61:2187-2194, 2012 




Type 1 diabetes has continuously increased world- 
vwde, and the highest incidence is found in Finland 
(1). Diabetic neptiropathy (DN) is a complication 
that develops in approximately 25^0% of patients 
with type 1 and type 2 diabetes. DN is the leading cause of 
end-stage renal disease (ESRD) in the developed world. 
Currently, 44% of the new cases of ESRD in the U.S. annually 
are attributable to DN (2). A better understanding of the 
causal factors of DN and its pathogenesis may lead to new 
strategies to decrease its incidence, preemptively treat the 
disorder, attenuate morbidity and mortality, and would be 
a valuable contribution to global public health. 

Several observations suggest that DN, one of the mgjor 
complications of type 1 and type 2 diabetes, has an in- 
herent genetic susceptibility. Familial clustering of DN is 
evident for both type 1 and type 2 diabetes (3-6), and ge- 
netic risk factors are being sought in multiple populations 
(7-9). Unfortunately, robust replication of many initial 
associations has not been forthcoming (10). 

This study recruited a large collection of individuals with 
lype 1 diabetes as part of the GEnetics of Nephropathy-an 
International Effort (GENIE) consortium and examined se- 
lected candidate loci associated with DN from genome-wide 
case-control studies or other association studies that 
reported high levels of statistical significance. The variants 
examined and the rationale for their inclusion are as follows: 

i) A single nucleotide polymorphism (SNP) (rsl617640) 
within the promoter region of the EPO gene (encoding 
erythropoietin) was identified as having a genome-wide 
significant (P < 5 X 10~^) association with ESRD and 
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proliferative diabetic retir\opathy (PDR) (11). Interest- 
ingly, erythropoietin levels were elevated sevenfold 
in the human vitreous fluid of nondiabetic individuals 
with the risk genotype TT compared with those with 
the wild-type GG genotype. In addition, EPO expres- 
sion levels were significantly elevated above control in 
the tissues and vitreous fluid of animal models of DN 
(DN in db/db mice) and in proliferative retinopathy 
(murine oxygen-induced retinopathy model), respec- 
tively (11). 

2) The engulfment and cell motility 1 gene (ELMOl ) has 
been reported to be associated with DN in Japanese 
patients with type 2 diabetes (12). Recently, Pezzolesi 
et al. (13), using the Genetics of Kidneys in Diabetes 
U.S. Study (U.S. GoKinD) cohorts, also examined 
ELMOl for association with DN and presented evi- 
dence of association of variants within this gene for 
the development of DN. However, the risk alleles for 
ELMOl identifled in their study differed from those 
reported in the original Japanese investigation. In the 
context of a genome-wide association study (GWAS), 
118 SNPs were assessed in 1,705 individuals of Euro- 
pean ancestry with type 1 diabetes (885 control subjects 
and 820 DN case subjects). The strongest associations 
in ELMOl in the U.S. study occurred at rsll769038 
(odds ratio [OR] 1.24; P = 1.7 X 10"^) and rsl882080 
(OR 1.23; P = 3.2 X 10"^), located in intron 16. Two 
additional SNPs, located in introns 18 and 20, were 
also nominally associated with DN. In total, eight 
ELMOl SNPs were reported to confer risk for DN, al- 
though none reached genome-wide signiflcance (13). 
Supportive evidence was also found in African Ameri- 
cans with type 2 diabetes and ESRD (14). 

3) The U.S. GoKinD GWAS analyzed 359,193 SNPs in 820 
case subjects (284 with proteinuria and 536 with ESRD) 
and 885 control subjects with type 1 diabetes but no 
evidence of DN. Although no risk variant achieved 
genome-wide signiflcance, the primary association anal- 
ysis identifled 11 SNPs representing four distinct chromo- 
somal regions (P < 1 X 10~^). The strongest association 
with DN reported in this study was on chromosome 9q 
with IS10868025 (OR 1.45, P = 5.0 X 10"^) (15). 

4) Finally, in an effort to systematically explore and compre- 
hensively capture common genetic variations that might 
be associated with DN, we reviewed the largest meta- 
analysis published to date studying genetic associations 
with the DN phenotype (7). In GENIE, we examined the 



top-reported SNP (or proxy) for each gene in that report 
for an association with DN. 

In this study, we have assembled the largest reported 
case-control sample of DN in type 1 diabetes to evaluate 
the previously reported genetic associations in newly geno- 
typed samples from the U.K., Republic of freland (R.O.I.), 
and Finland, plus pre-existing data from the U.S. GoKinD. 



RESEARCH DESIGN AND METHODS 
Cohorts 

U.K.-E.O.I. collection. Recruited individuals were part of the All Ireland- 
Warren 3-Genetics of Kidneys in Diabetes U.K. collection (U.K.-R.O.I.). All were 
self-reported as white, with grandparents bom in the U.K. or beland, and type 1 
diabetes diagnosed before the age of 31 years requiring uninterrupted insulin 
treatment. Case subjects (n = 903) with DN had persistent proteinuria (>0.5 
g/24 h), hypertension (> 135/85 mmHg and/or treatment with antihypertensive 
medication), and diabetic retinopathy. ESRD (28%) was defined as requiring 
renal replacement therapy or having received a kidney transplant hidividuals 
in the control group (n - 1,001) had had type 1 diabetes for at least 15 years, 
had no evidence of microalbuminuria on repeated testing, and were not re- 
ceiving antihypertensive medication (Table 1). 

Finnish Diabetic Nephropathy Study (FinnDiane ). The FinnDiane study 
is a nationwide multicenter study of >4,800 adult participants with type 1 
diabetes (16). This study comprises genotype data for 2,914 patients with type 
1 diabetes diagnosed before age 35 years and insuhn treatment started within 
1 year of diagnosis. The disease status was defined by urine albumin excretion 
rate (AER) or urine albumin-to-creatinine ratio (ACR) in at least two of three 
consecutive urine collections at local centers. Macroalbuminuria (n - 686) 
was defined as AER >200 (jLg/min or >300 mg/24 h or an ACR >25 mg/mmol 
for men and >35 mg/mmol for women in overnight, 24-h, or spot urine col- 
lections, respectively. Similarly, the limit for normal AER (n = 1,601) was <20 
(ig/min or <30 mg/24 h or ACR <2.5 mg/mmol for men and <3.5 mg/mmol for 
women. Control patients with normal AER were required to have type 1 
diabetes duration of at least 15 years. ESRD (n = 627) was defined as ongoing 
dialysis treatment or receipt of a kidney transplant. From the total, 505 par- 
ticipants were included from an independent Firmish cohort collected by the 
National Institute for Health and Welfare (17). These participants met the 
FinnDiane diagnosis and selection criteria and were analyzed together with 
the FiimDiane cohort (Table 1). 

U.S. GoKinD. The U.S. GoKinD study consists of a DN case-control cohort of 
individuals diagnosed with type 1 diabetes before age 31 years, who were 
between 18 and 59 years of age at enrollment, and who began insulin treatment 
within 1 year after diagnosis (18). The 905 case subjects were defined as people 
aged 18-54, with type 1 diabetes for at least 10 years, and DN. The 898 control 
subjects were aged 18-59, had type 1 diabetes for at least 15 years, but did not 
have DN. The DN definition includes individuals with ESRD (on dialysis or 
having received a kidney transplant) or persistent macroalbuminuria (at least 
two of three tests positive for albuminuria by dipstick >1+ or ACR >300 ^.g 
albumin/mg of urine creatinine). The U.K. GoKinD inclusion criteria were used 
to recruit individuals to the control group, hidividuals were recruited at two 
study centers, George Washington University (GWU) and the Joslin Diabetes 



TABLE 1 

Phenotypic characteristics of the GENIE cohorts (U.K.-R.O.I., FinnDiane, and U.S. GoKinD) 





U.K.-R.O.I. 


FinnDiane 


U.S. GoKinDt 




Case subjects 


Control subjects 


Case subjects 


Control subjects 


Case subjects 


Control subjects 


n 


903 


1,001 


1,289 


1,577 


774 


821 
















Male 


531 


438 


764 


651 


402 


342 




Type 1 diabetes 


Duration (years) 


32.9 ± 9.5 


27.2 ± 8.7 


32.8 ± 9.1 


27.9 ± 9.5 


31.4 ± 7.8 


25.4 ±7.7 


Age at diagnosis (years) 


14.6 ± 7.7 


14.5 ± 7.8 


12.8 ± 7.6 


15.1 ± 8.3 


11 ± 6.6 


13 ± 7.3 


HbAie (%) ■■■i 


1 9.0 ± 1.9 


8.7 ± 1.6 


8.8 ± 1.6 


8.0 ± 1.2 


7.5 ± 1.9 


7.5 ± 1.2 


BMI (kg/m^) 


26.2 ± 4.7 


26.2 ± 4.2 


25.5 ± 4.2 


25.2 ± 3.4 


25.7 ± 5.2 


26.1 ± 4.3 


ESRD (%) 


28.0 


0 


48.2 


0 


65.6 


0 



Categoric data are shown as indicated; continuous data as mean ± SD. fReanalysis of the U.S. GoKinD dataset using new quality control 
filters to account for published plate effects (see resbasch design and methods for complete details). 
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Center (JDC) using differing methods of ascertainment and recruitment (see 
Pezzolesi et al. [15] for details). Analysis of the U.S. GoKinD cohort was lim- 
ited to individuals whose reported primary ethnicity was white. 
Phenotype definition: DN and EPO study outcomes. DN was the primary 
outcome for all association studies oftheSNPs investigated. For the EPO study, 
we used the phenotypic definitions used in the original report (11); namely, 
ESRD as defined by case subjects who were on dialysis or who had received 
a renal transplant; concurrently, this ESRD population also had to have evi- 
dence of advanced PDR on physical examination and a history of laser 
treatment. Control subjects with evidence of PDR were excluded. We also 
examined PDR as an outcome independent of ESRD status. For this analysis, 
case subjects had clinical evidence of PDR, whereas control subjects had none 
(irrespective of DN status). Cohorts analyzed in the original EPO study (Fig. 1) 
were composed of European-American cases, and control subjects were col- 
lected from distinct geographic areas in the U.S. These included the GoKinD 
cohort (Boston, Pittsburgh, and Minnesota), the Utah cohort (Salt Lake City), 
and the Boston cohort (Boston Joslin Center for Diabetes) (11). 
SNP selection. SNP markers with evidence for association with DN sus- 
ceptibility in reference studies (11) were selected for genotyping in GENIE. 
Where more than one SNP was associated at a particular locus with DN, the 
most strongly associated variant was selected for genotyping in the DN 
case-control cohorts. Where no genotyping assay could be developed for the 
index SNP, a proxy in strong linkage disequilibrium (LD) was genotyped using 
the CEU HapMap population. CEU is the official three letter code for the 
HapMap samples of Utah residents with ancestry from Northern and Western 
Europe (see http://hapmap.ncbi.nlm.nih.gov/citinghapmap.htral). For the 
original U.S. GoKinD reported results and ELMOl, additional SNPs within 
20-kb upstream and downstream of the locus (or index SNP) were selected 
using the SNP Annotation and Proxy Search (http://www.broadinstitute.org/ 
mpg/snap/), specifying chromosome position, CEU samples, and 1000 Genome 
Pilot 1 data. The expanded SNP list was extracted from GWAS results for 
U.K.-R.O.I. and FinnDiane (N. Sandhohn et al., submitted). 
De novo genotyping. For the U.K-R.O.I. collection (n = 1,904 unique indi- 
viduals), SNPs were genotyped using Sequenom iPLEX (Sequenom, Hambxirg, 
Germany) or TaqMan (AppUed Biosystems, Warrington, U.K.) technology. 
Duplicate and no-DNA-template samples were included on all plates as ex- 
perimental controls. 

In FiimDiane, the EPO locus was genotyped with TaqMan chemistry (Applied 
Biosystems, Foster City, CA) in 3,363 samples, of which 251 case subjects 
(ESRD + laser treatment) and 987 control subjects (no DN, no retinopathy) 
passed the phenotype criteria. All other SNPs were genotyped at the Institute of 
Molecular Medicine Finland (Helsinki, Finland) on the Dlumina's BeadArray 
OlOQuad array. Dlumina's BeadStudio clustering algorithm was used to call 
genotypes in FinnDiane. SNPs were filtered for those with call rates >95%, minor 
allele frequency (MAF^ >1%, and test for Hardy Weinberg equilibrium (HWE; 
P > 1 X 10"^). Sample filters included individual call rates >95% and no first- 
degree relatives and resulted in 1,289 case subjects and 1,577 control subjects. 

Existing genotype data for the U.S. GoKinD genotype data were dovpnloaded 
from dbGAP (phs000018.v2.pl, retrieved June 2010), containing genotype data 
from the Affymetrix 500 K array (Affymetrix, Santa Clara, CA). The version 2 
genotype data differed from the original U.S. GoKinD data, containing updated 
and recalled genotype calls for a previously reported problematic plate (19), 
and SNP call rate (>90%), MAF (>1%), and HWE performed by the National 
Heart, Lung, and Blood Institute. We applied additional qualify control (QC) 
filters to the data, including removing samples with evidence of contamination 
(extreme heterozygosify, n = 16), known parents (n = 4), non-European- 
reported ancestry (n = 121), reported sex-genotype mismatch (w = 1), cryptic 
related subjects (n = 4), and principal component analysis admixture outliers 
(n = 5). Furthermore, SNP QC filters were appUed to remove known prob- 
lematic SNPs, HWE (P > 1 X 10"''), call rate (90%), missingness, MAF (1%), 
and SNPs plate effects. 

GWAS genotyping. Genotypes for U.K.-R.O.I. (n = 1,830) and FinnDiane 
(w = 3,651) were supplemented with SNPs retrieved from GWAS results for 
each study (Sandholm et al., submitted). The FiimDiane study samples were 
genotyped using the Dlumina's BeadArray 610-Quad (Illumina, San Diego, CA) 
array at the Institute of Molecular Medicine Finland (FIMM, Helsinki, Finland). 
The U.K-R.O.I. study samples were genotyped at the Broad Institute using the 
Dlumina Omnil-Quad array. Samples with insufficient DNA quahfy, quantify, 
or poor genotype concordance with previous available genotypes at a finger- 
printing stage were excluded. For all three discovery datasets (FinnDiane, 
U.K.-R.O.I., U.S. GoKinD), a standardized and detailed genotype QC proce- 
dures were applied using PLINK supplemented with perl and R scripts. These 
included selecting SNPs with a call rate >90%, MAF >1%, and concordance 
with HWE (P > 10"^). Wc discarded samples with a call rate <95%, extreme 
heterozygosify, cryptic relatedness, or ethnic outliers by principal components 
analysis, and tested for SNP missingness by haplotype (P < 10~^, by phenotype 
(P > 10"^), and by plate effects (P < 10 '). Probes for copy number variation 



as weU as sex chromosome and mitochondrial SNPs were excluded from 
analyses. After applying the QC protocol, we had access to 549,530 SNPs in 
3,370 FinnDiane samples, 791,687 SNPs in 1,726 U.K-R.O.L samples, and 
360,899 SNPs in 1,595 U.S. GoKinD samples. 

SNP imputation. MACH 1.0 software (ht^:/Avww.sph.uiracKedu/csg/abecasis/ 
MACH) with the HapMap phase n CEU reference panel was used to perform 
SNP imputation for GWAS results in each cohort. Estimates of the crossover 
and error rates were obtained via 50 iteration rounds in —300 randomly se- 
lected samples per cohort A greedy algorithm was used for imputation, and 
the maximum likelihood method was specified to yield allele dosages. A filter 
was applied to exclude SNPs with low imputation quahfy (r^ < 0.6), 
resulting in ~2.4 million SNPs per cohort. 

Statistical analysis. Association tests were conducted using PLINK vl.07 (20) 
(http://pngu.mgh.harvard.edu/purcell/plink), with logistic regression acjjusted 
for sex and age. U.K.-R.O.L was a^usted for recruitment center, but the two U.S. 
GoKinD centers, GWU and JDC, were analyzed separately as reported by 
Pezzolesi et al. (15). Data from the GWAS genotyping was adjusted addition- 
ally for duration of type 1 diabetes and principal components from Eigenstrat 
analysis. The EPO locus was analyzed with the Pearson test in the Finn- 
Diane dataset without adjusting for any covariates. Fixed-effects meta-analyses 
were conducted with the software package Comprehensive Meta-Analysis 
(Version 2.2.040, Englewood, NJ) and the software package METAL (http:// 
www.sph.umich.edu/csg/abecasis/Metal/) (21) under the additive genetic 
model To determine the appropriate statistical cutoff for correction for 
multiple testing, we calculated the total number of effective tests (because 
a large portion of SNPs were in LD), using SNPSpD (http://gump.qimr.edu.au/ 
general/daleN/SNPSpD/). SNPSpD uses correlation between analyzed SNPs to 
calculate the total number of independent tests (22). The total number of SNPs 
tested is 2,199, with 113.7 effect-independent tests. Thus the experiment-wide 
cutoff for statistical significance was set at 4.4 x 10~* (0.05/113.7). 



RESULTS 

EPO promoter polymorphism. The association of the 
EPO promoter polymorphism, rs 1617640, with DN was 

evaluated by de novo genotyping of the SNP in GENIE. For 
this analysis, case subjects were defined as having both 
ESRD and PDR because the initial report showed the 
polymorphism was robustly associated with DN when both 
of these "extreme" phenotypes were coexpressed. Significant 
association was not observed in the U.K.-R.O.I. (P = 0.19) 
or FinnDiane collections (P - 0.60), although the direc- 
tions of effect were consistent with the original report. 
Fixed-effects meta-analysis of the association of rsl617640 
with ESRD/PDR, including the previously reported cohorts 
(a total of 3,162 case and 3,845 control subjects across five 
separate cohorts of European and European-American 
ancestry) retained genome-wide statistical significance 
(OR 1.31 [95% CI 1.20-1.44], P = 2 X 10"''; Fig. 1). 

As an additional experimental control, we examined the 
potential association of the EPO promoter polymorphism 
with the development of PDR in case subjects, irrespective 
of ESRD status. No association was observed between 
EPO and PDR for the individual cohorts or in the meta- 
analysis of the combined results for FinnDiane (OR 0.95 
[95% CI 0.85-1.04], P = 0.25) or and U.K.-R.O.L (0.96 [0.88- 
1.04], P = 0.29). Furthermore, no association was observed 
after inclusion of the restricted phenotype, PDR, in U.S. 
GoKinD case and control subjects, separately, or in the 
meta-analysis of all cohorts combined (results not shown). 
ELMOl . Neither single-center nor meta-analysis of de novo 
genotyping in U.K.-R.O.L, nor GWAS data for FinnDiane, 
revealed a significant association in subjects with type 1 
diabetes between rs741301, the previously reported risk 
variant within ELMOl, and DN (OR 1.04 [95% CI 0.95-1.13], 
P = 0.46; Supplementary Fig. 1). Pezzolesi et al. (13) also 
tested rs741301 but did not replicate the reported associa- 
tion. They went on to test other SNPs in the region and 
reported nominal associations (P = 0.002-0.05) with eight 
other SNPs . We examined LD between rs741301 and the 
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EPO rs1617640 and Risk of Combined PDIR/ESRD 



study name SNP 



Statistics for each study 



Odds Lower Upper 

ratio limit limit 

Utah rs1 6 17640 1.446 1.146 1.825 

US GoKinDrsI 61 7640 1.535 1.320 1.785 

Boston rs1617640 1.382 1.050 1.819 

Finn Diane rs161 7640 1.055 0.866 1.285 

UKROI rs161 7640 1.140 0.938 1.386 

1.313 1.202 1.435 



Z-Value p-Value 

3.107 0.002 



5.566 
2.308 

0.532 
1.315 
6.027 



0.000 
0.021 
0.595 
0.188 
2x10-8 



Odds ratio and 95% CI 



0.5 



o 



Minor Allele Major Allele 



Meta Analysis 

FIG. 1. Previously published and new results from this study provide an estimate of the effect of the EPO promoter SNP (rsl617640) on the risk of the 
combined phenotype of FDR and ESRD in type 1 diabetes in five cohorts (3,162 case subjects and 3,845 control subjects) in a fixed-effects meta-analysis. 



SNPs reported to be associated with DN by these inves- 
tigators. The statistic betweer\ rs741301 and the other 
SNPs revealed only low to moderate LD, ranging from 0.38 
to 0.65 (Supplementary Table 1). 

As an additional and more extensive test of variants in 
this region, we performed an expanded analysis capturing 
all available SNPs 20 kb upstream and downstream of the 
ELMOl locus to account for LD differences presumably 
due to ancestry. Results of the expanded analysis did not 
reveal any significant SNPs for either cohort individually 
or for the two cohorts meta-analvzed after correcting for 
multiple testing (P < 4.3 X Furthermore, no SNP 

achieved significance with inclusion of the U.S. GoKinD 
results in the meta-analysis (Supplementary Table 2). 
Rislc variants reported from U.S. GoKlnD for DN in 
type 1 diabetes. Eleven DN susceptibility SNPs that were 
first reported by the U.S. GoKinD investigators as highly 
associated with the risk of developing DN in type 1 di- 
abetes were parsed into eight candidate loci. We selected 
one representative SNP for each region of strong LD in 
which multiple SNPs represented the same association 
signal. After performing additional QC checks (see research 
DESIGN AND METHODS), we first tested the eight U.S. GoKinD 
potential DN susceptibility SNPs by reanalyzing the U.S. 
GoKinD dataset downloaded from dbGAP (19). As in the 
original report, none of the eight SNPs were associated 
with DN at genome-wide statistical significance (Table 2). 
The two SNPs in the FERM (F [Band 4.1], E [Ezrin], R 
[Radixin], M [Moesin]) domain 3 (FRMD3) region showed 
similar P values as in the original report (2.1 X 10~^ and 
1.6 X 10~^, respectively; Table 2). In our analysis, the 
statistical significance of the P values for SNPs in the 
CPVL/CHN2 and CARS regions was reduced from 6.5 X 
lO"'^ to 2 X 10"^ and 6.4 X 10"^ to 2.2 X 10"^, respectively. 
P values for the 6 SNPs in the 13q region were also 
changed from 1.8-7.0 X 10"*^ to 1.4-9.5 X 10"^. 

We next examined these SNPs in our newly genotyped 
samples. Case-control association analysis for these loci in 
the U.K.-R.O.I. and FinnDiane samples revealed no signifi- 
cant associations in either cohort. The strongest signal was 
observed at rs39075 near CPVL/CHN2 in U.K.-R.O.L, (OR 
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1.12, P = 0.08) and in meta-analysis of the two replication 
cohorts (U.K.-R.O.I. and FinnDiane; OR 1.06, P = 0.06; 
Table 2). In expanded locus-region analyses (plus or minus 
20 kb of the locus of interest), no SNP reached significance 
after ac^ustment for multiple testing (one-tailed P = 0.03, 
experiment-wide threshold P = 4.3 X 10~*) for the two 
cohorts separately or via meta-analysis. The combined 
meta-analysis including U.S. GoKinD revealed two SNPs 
downstream of FRMD3, rsl888747 {P = 1.5 X 10"^) and 
rsl3288659 (P = 9.7 X 10"^), which showed significance 
after acljusting for experiment-wide multiple testing (P < 
4.3 X 10~*). However, neither SNP achieved genome-wide 
significance (P < 5 X 10~^; Supplementary Table 3). 
Pooled meta-analyses examining variants associated 
with DN in type 1 and type 2 diabetes. In the most 
comprehensive literature search for DN associated genetic 
variants to date, Mooyaart et al. (7) identified 24 loci. In 
GENIE, we examined all the available top-reported SNPs 
(or their proxies) for each gene in that report. Three SNPs 
were nominally associated (P < 0.05) with DN: rsl3293564 
at UNC13B (P = 0.01) and rsl79975 at the ACE (jP = 0.03) 
in FinnDiane, and rs39075 at CPVL/CHN2 (P = 0.05) in the 
U.K.-R.O.L samples. In a meta-analysis of the two cohorts, 
the ACE polymorphism remained nominally significant 
(P = 0.04). Including the U.S. GoKinD results, the FRMD3 
signal at rs 1888747 emerged as noted above; no other 
signals were significant after adjusting for multiple com- 
parisons (see Supplementary Table 4 for ftdl details). 



DISCUSSION 

Using a large, homogeneous population sample of European 
ancestry subjects with type 1 diabetes in the GENIE con- 
sortium, we were imable to replicate most of the previously 
reported genetic associations with DN that we examined. 

Our findings do not support previously reported genetic 
associations with DN in type 1 diabetes in the largest GWAS 
published to date (15). None of those signals reached 
genome-wide statistical significance with the addition of 
larger, similarly ascertained datasets. Using ORs at the 
lower limit of the 95% CI from the original publication, our 
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combined sample of U.K. -R.O.I, and FinnDiane cohorts in 
GENIE had 99.9% power (at a = 0.001) to detect the U.S. 
GoKinD reported effect sizes (Supplementary Table 5); 
thus, our negative results (even at a = 0.05) were un- 
expected. Further, after performing additional QC checks 
on the original U.S. GoKinD dataset and combining these 
samples with GENIE, we were unable to achieve genome- 
wide significant replication in the meta-analysis. 

We investigated other previously reported genetic asso- 
ciations with DN in subjects with type 1 diabetes, namely 
the EPO promoter polymorphism rsl617640 and variants 
within the ELMOl gene; these also failed to replicate in 
GENIE. The EPO promoter polymorphism retained genome- 
wide significance after meta-analysis of the prior data 
combined with ours, for the combined phenotype of ESRD 
and PDR. However, the overall P value attained was at- 
tenuated from 2.8 X 10"" to 2 X 10"*'. 

We also did not observe evidence for replication of the 
association of rs741301 in the ELMOl gene and DN in 
GENIE. Shimazaki et al. (12) first reported an association 
for this genetic variant with DN in Japanese subjects with 
type 2 diabetes. Subsequently, Pezzolesi et al. (13) failed to 
rephcate the association of the same SNP in type 1 diabetes 
but reported that eight SNPs within the gene locus had 
nominal associations with DN in the European-derived 
subjects in U.S. GoKinD. Lack of replication may be due to 
different genomic patterns in populations of diverse ances- 
tries or the distinct genetic architecture of DN in type 1 
versus type 2 diabetes (14). To address the first concern, 
we expanded our analysis by interrogating all available 
SNPs plus or minus 20 kb of ELMOl. Because in the U.S. 
GoKinD analysis of ELMOl the risk variants were only in 
weak to modest LD with the index SNP (r^ between 
rs741301 and the U.S. GoKinD SNPs ranged from 0.38 to 
0.65), we reasoned that this expanded strategy would ac- 
count for most of the differences in LD presumably due to 
ancestry. The results of this expanded analysis, however, 
did not reveal any significant association for the additional 
670 SNPs interrogated in ELMOl in the U.K.-R.O.I. or 
FinnDiane cohorts individually, for these two cohorts in 
meta-analysis, or when combined with previously reported 
risk variants reported from U.S. GoKinD. 

As a final step to ensure a systematic and more com- 
prehensive approach to candidate loci associated with DN, 
we reviewed the top SNPs from the 24 loci cited by the 
largest meta-analysis published to date that has examined 
candidate genetic variants for association with DN (7). From 
the results of the pooled meta-analyses, the strongest signal 
emerged at rsl888747 in FRMD3 (P = 1.5 X 10"*). Variants 
within this gene have been reported to be associated with 
DN in European cohorts (7,15,23), and Freedman et al. 
(24) recently reported evidence for a role for gene-gene 
interactions between myosin heavy-chain 9 (MYH9)- 
apolipoprotein LI (APOLl ) haplotypes and FRMD3 in 
African American subjects. It should be noted, however, 
that risk variants within this gene have achieved only 
nominal but not genome-wide statistical significance in all 
previous reports. 

From the foregoing observations, we conclude that 
there is a high Ukelihood that many of the previously 
reported positive associations with DN represent potential 
false-positive findings (type I error). We emphasize that the 
combined sample of U.K.-R.O.I. and FinnDiane represents 
a substantially larger collection of case and control subjects 
than U.S. GoKinD and is well powered to detect the origi- 
nally reported effect sizes, thereby decreasing the likelihood 
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of a false-negative finding (type II error), even accounting 
for the li]<ely overestimation of effect sizes due to the 
winner's curse phenomenon (25). We also harmonized the 
ascertainment criteria for case-control definitions across 
all the study populations (including U.S. GoKinD), making 
it unlikely that phenotypic heterogeneity across study 
populations explains the lack of replication. 

A crucial issue that bears on the interpretation of 
case-control studies of the genetics of DN concerns the 
adequacy of phenotype definition. In this and most studies 
cited to date, there is the presumption that long-duration 
diabetes exposure and the presence of frank protein in the 
urine — macroalbuminuria — defines DN and that pheno- 
typic heterogeneity has been well controlled through this 
classification. These definitions are derived in large mea- 
sure from the classic studies of Parving et al. (26), Viberti 
et al. (27), and Mogensen and Christensen (28), who docu- 
mented 30 years ago a virtually inexorable progression to 
ESRD in patients who developed microalbuminuria after 
approximately 2 decades of exposure to the diabetic met- 
abolic milieu. However, these longitudinal findings were 
based on small numbers of patients. The plasticity of DN 
phenotypes is reflected in more recent and much larger 
longitudinal studies showing that most patients with type 1 
diabetes, categorized initially as having microalbuminuria, 
undergo regression to normoalbuminuria with preservation 
of renal function (29). It is not entirely clear that micro- 
albuminuria versus macroalbuminuria, stage of chronic 
kidney disease and attendant renal function, the rate of 
renal decUne, or the occurrence of extreme phenotypes, 
such as ESRD/PDR, represent one disease process along 
a continuum or many distinct disease states, each of which 
may be under distinct genetic control. As pointed out re- 
cently (24), genetic variants, such as those in MYH9 and 
APOLl that axe common in certain ethnic groups, may 
mask the effects at other loci unless methods such as 
muItUocus modeling and interaction analyses are used to 
control for these effects. 

In addition, phenotypic variation may be a function of 
ethnicity and disease-specific gene expression. For exam- 
ple, Pima Indians with type 2 diabetes have very early- 
onset DN, characterized by an accelerated loss of renal 
function and progression to ESRD despite lower blood 
pressures and lipid levels, factors thought to be protective 
(30). This has been postulated to be due to structural dif- 
ferences in the nephron-podocyte number and density per 
glomerulus ("podocyte insufficiency"), a decrease in net 
nephron mass (glomerulopenia) resulting in glomerulomegaly, 
increased intraglomerular capillary pressure, and ultimately, 
hyperfiltration injury (30). Whether these structural and 
intrarenal hydratilic changes could be genetically regulated 
is ultimately a testable hypothesis; they warrant further in- 
vestigation to continue the inquiry why certain populations 
have an apparent disproportional susceptibility to ESRD 
and, particularly, DN. 

In summary, we have presented evidence that several 
previously reported genetic associations with DN in type 1 
diabetes could not be replicated in a large, homogeneous 
sample of subjects with type 1 diabetes. Oiu: failure to 
replicate these associations underscores the need to apply 
stringent statistical thresholds of significance, maximize 
power through meta-analysis of all available data, and seek 
replication in independent samples, as has been proposed 
by a number of different authors (31,32). Finally, the ap- 
plicability and generalizability of DN risk loci from type 1 
diabetes to type 2 diabetes, and the related question of 



shared genetic susceptibility for nephropathy between 
type 1 and type 2 diabetes, remain imresolved. 
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